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Abstract. In this paper we consider the geodesic X-ray transform with attenua- 
tion coefficient as a combination of smooth complex function and 1-form. We show 
that attenuated X-ray transform applied to the pair of tensors is injective modulo 
the natural obstruction. 



1. Introduction 

The geodesic X-ray transform of a function is denned by integrating over geodesies. 
It is naturally arises in linearization of the problem of determining a conformal factor 
of a Riemannian metric in the boundary rigidity problem. The X-ray transform also 
arises in Computer Tomography, Positron Emission Tomography, geophysical imaging 
in determining the inner structure of the Earth, ultrasound imaging. Uniqueness 
result and stability estimates of the geodesic X-ray transform were obtained by R. 
G. Mukhometov [16] for simple surface. For simple manifolds of any dimension this 
result was proven in [HE], see also V. A. Sharafutdinov's books [25] [26] . In his paper 
N. S. Dairbekov generalized this result for nontrapping manifolds without conjugate 
points [6]. Fredholm type inversion formulas were given in [22] by L. Pestov and G. 
Uhlmann. 

In our paper we consider the uniqueness problem for the attenuated X-ray transform 
on a surface. Now let us define this transform. 

Let M be a surface with smooth boundary dM and endowed with Riemannian 
metric g. 

Define the sets of inward and outward unit vectors respectively by 
d±SM = {(x, £) e SM : x e dM, ±(f , v{x)) > 0}, 

here and futher v is the unit inner normal to dM. 

We assume that (M, g) is nontrapping, which means that the time t(x,£) when the 
geodesic y x ,t(t) exits is finite for each (x,£) € d + SM. 

Let h G C°°(M) and a be a smooth 1-form on M. Consider an attenuation coeffi- 
cient o as a combination of h and a, i.e. a(x, £) = h(x) + at x (Q f° r i x i£) e SM. Let 
ip : SM — > R be a smooth function on SM. Define the attenuated X-ray transform 
of ip by 



I a ip(x, := / ^(7(t),7(t))exp / a(y x>i (s), 7 x>€ (s)) ds 

"T(*,£) 



t 

' dt. 



o 







t 

' dt. 



^(7(*),7(*))exp / h(j x ^(s)) + a Jx ^ s )(%^(s)) ds 
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where (x,g) G d + SM. 

In this paper we will consider the problem of injectivity of attenuated X-ray trans- 
form applied to functions on M. 

We say that (M, g) is simple if 1) DM is strictly convex and 2) all geodesies have no 
conjugate points. In this case, M is diffeomorphic to the unit ball of IR 2 . In particular, 
a simple manifold is nontrapping. The notion of simplicity arises naturally in the 
context of the boundary rigidity problem [15j . 

The Euclidean attenuated X-ray transform is the basis of the medical imaging 
modality SPECT. The attenuated geodesic ray transform arises in inverse transport 
problems with attenuation [T51 [TJ] , when the index of refraction is anisotropic and 
represented by a Riemannian metric. It also arises in geophysics where there is 
attenuation of the elastic waves. Rather unexpectedly, this transform also appeared 
in the recent works [H [T_2] in the context of Calderon's inverse conductivity problem 
in anisotropic media. 

The following result shows that attenuated X-ray transform on 2-dimensional sim- 
ple Riemannian manifold is injective. 

Theorem A. Assume (M, g) is a simple Riemannian surface. Let h be a smooth 
complex function and let a be a smooth complex 1-form on M . Consider an atten- 
uation coefficient a as a combination of h and a, i.e. a = h + a. Suppose that 
f : M — > C is smooth such that I a f = ; then f = 0. 

In the case where M = M 2 with the Euclidean metric, the corresponding injectivity 
result for the attenuated X-ray transform with a = has been proved by different 
methods in E. V. Arbuzov, A. L. Bukhgeim and S. G. Kazantsev [3], R. G. Novikov 
[19], F. Natterer |18j . and J. Boman and J.-O. Stromberg [5]. These methods also 
come with inversion formulas. If M is the unit disc in IR 2 with Euclidean metric, a 
direct inversion formula was given by Kazantsev and A. A. Bukhgeim [11] for the 
case a = 0. See D. V. Finch |8J and P. Kuchment [10] for surveys of these and 
other developments in Euclidean space. Injectivity of the attenuated X-ray transform 
for simple surface has been proved in |24J by M. Salo and G. Uhlmann when a is a 
function on M and in [2U] by G. P. Paternain, M. Salo, G. Uhlmann when a is purely 
imaginary. 

Now we are interested in I a applied to the functions on SM of the following type 

Vtofl = f, ; ...,J.rK l! ■■■€- + ft.-^Wf 1 • ■■C m -\ 

where / and /3 are smooth symmetric (covariant) tensors of ranks m and m — 1 on 
M. It has nontrivial kernel since 

I a (hp + d s p + a(pa)) = 

for any smooth symetric (m — l)-tensor p with p\qm = 0. Here and futher a denotes 
symmetrization. We investigate if these are the only elements of the kernel. 

The case m — 1 corresponds to the geodesic weighted Doppler transform in which 
one integrates a vector field along geodesies. This transform appears in ultrasound 
tomography to detect tumors using blood flow measurements and also in non-invasive 
industrial measurements for reconstructing the velocity of a moving fluid. For m > 1, 
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problems of this kind also appear in inverse problems for hyperbolic equations [23] . 
photoelasticity pQ. The integration of tensors of order two along geodesies when a = 
arises as the linearization of the boundary rigidity problem and the linear problem is 
known as deformation boundary rigidity. 

The next theorem describes elements in the kernel of the attenuated X-ray trans- 
form applied to combination of smooth symmetric tensors of ranks m and (m — 1). 



Theorem B. Let (M,g) be a simple 2- dimensional Riemannian manifold, and con- 
sider h and a to be a smooth complex function and 1-form (resp.) on M . Denote 
a = h + a. If f and (3 are smooth symmetric m and (m — l)-tensors (resp.) on M 
such that I a (f + (3) = 0, then f = hp and (3 = d s p + a (pa) for some smooth symetric 
m — 1-tensor p with p\qm = 0. 

For m = 1 results of this type were given in the unit disc in IR 2 in [TT] when a = 0, 
for simple manifolds with ||a||jr,oo small and a = in [7], and for simple manifolds 
with g and a close to real analytic in [9J. Corresponding result have been proved in 
[20] by G. P. Paternain, M. Salo, G. Uhlmann when a is purely imaginary. 

When attenuation coefficient vanishes this is the case of the ordinary X-ray trans- 
form. Some earlier results corresponding to this problem were described in [21J and 
Theorem B was proved for the case when attenuation coefficient is absent. 

Proof of Theorem A follows the same scheme as in [21] and proof of Theorem B 
is mimic those in J2JJ. Briefly, the paper is organized as follows. In Section 2 we 
recall some facts and definitions from [20, 24J that will be used in our paper. These 
are notions of (anti)holomorphic functions, existence of (anti)holomorphic integrating 
factors for the transport equation. Next we recall regularity of solutions of transport 
equation. In Section 3.1 we prove injectivity of the attenuated X-ray transform using 
pseudodifferential arguments. Finally, Section 3.2 is devoted for the proof of Theorem 
A. Here we use the same method as in the second proof of Theorem 1.1 in [21] see 
Section 5]. 



2. Preliminaries 

For (x,£) G SM we denote by : [t-(x,£),t+(x,£)] — > M a maximal geodesic 
with initial conditions j x ,z(Q) = x, %z(0) = £ and j x ^(t-(x,^)),j x ^(t + (x,^)) G DM. 
Using the implicit function theorem it is easy to show that the functions r_ and r + are 
smooth near a point (x, £) such that the geodesic j x &{t) meets DM transversally at 
t = T-(x, £) and t = t+(x,£) respectively. Since, by strictly convexity this condition 
holds everywhere on SM\S(dM), the functions r_ and r + are smooth on SM\S(dM) 
[251 Section 4.1]. 

The generator G of the geodesic flow <f> t is the vector field on SM defined as 

G = |> t (x,0 . 
at t=o 

Since M is assumed oriented there is a circle action on the fibers of SM with 
infinitesimal generator V called the vertical vector field. In local coordinates it is 
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given by 

The space L 2 (SM) decomposes orthogonally as a direct sum 

L 2 (SM) = @H k 



where H k is the eigenspace of — iV corresponding to the eigenvalue k. A function 
u G L 2 (SM) has a Fourier series expansion 

oo 

u= 22 Mfc > 

k=— oo 

where u k G H k . Also ||m|| 2 = Xlll^fcll 2 ' where ||u|| 2 = (m,m) 1//2 . The even and odd 
parts of u with respect to velocity are given by 

u + := 22 u ^ u - := ^2 Uk ' 

k even A; odd 

Locally we can always choose isothermal coordinates (x l ,x 2 ) so that the metric can 
be written as ds 2 = e 2X ((dx 1 ) 2 + (da; 2 ) 2 ) where A is a smooth real-valued function of 
x = (x 1 ^ 2 ). This gives coordinates (x l ,x 2 ,ip) on SM where p is the angle between 
a unit vector v and d/dx 1 . In these coordinates we may write V = 8/ dip and 

u k (x,<p) = Q^j u(x,t)e~ kt dty** = u k (x)e ik *. 

A function u on SM is called holomorphic if u k = for all k < 0. Similarly, we say 
that a function u on SM is called antiholomorphic if u k = for all k > 0. 

The following result that was obtained in f2M Lemma 2.6] will be used in our paper: 

Lemma 2.1. The product of two (anti) holomorphic functions is (anti) holomorphic, 
and e w is (anti) holomorphic if w is (anti) holomorphic. 

Now we recall some results about the existence of holomorphic and antiholomorphic 
integrating factors for the equation 

Gu + au = -ip in SM, 

where we assume that a(x, £) = h(x) + a x (0- All functions in this section will be 
scalar and complex valued. 

By a (anti) holomorphic integrating factor we mean a function e~ w , where w G 
C°°(SM) is (anti)holomorphic in the angular variable, such that for all r G C°°(SM) 

e w G{e- w r) = Gr + ar in SM. 

This is equivalent with the equation 

Gw = -a in SM. 

It was proved in [20, Theorem 4.1] that holomorphic and antiholomorphic integrat- 
ing factors always exist. This result is stated below. 
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Theorem 2.2. Let (M,g) be a simple Riemannian surface and a £ C°°(SM). The 
following conditions are equivalent. 

(a) There exist a holomorphic w £ C°°(SM) and antiholomorphic w £ C°°(SM) 
such that Gw = Gw = —a. 

(b) <x(x,£) = h(x) + a x (£) where h is a smooth function on M and a is a 1-form. 

Using this theorem we see that the boundary value problem 

Gu + au = —ip in SM, u\qsm = 
is equivalent with the problem 

(1) G(e- W u) = -e~ w i>, -e- w u\ dSM = 0. 

Note. If ijj — f(x) then the right hand side of ([1]) equals to e~ w f and is holomorphic 
since w is holomorphic. The next result also was proved in [2U Proposition 5.1] and 
tells that in this case e~ w u is holomorphic too. 

Lemma 2.3. Let (M,g) be a simple Riemannian surface and let f be a function on 
SM (anti) holomorphic in the angular variable. Suppose that v £ C°°(SM) satisfies 

Gv = -f in SM, v\ dSM = 0. 
Then v is (anti) holomorphic in the angular variable, and Vq = 0. 

Before finishing this section we state some regularity result for the transport equa- 
tion from [201 Proposition 5.2]. 

Lemma 2.4. Let F : SM -> C be smooth and I a (F) = then 

u{x,0= / F(0t(x,O)exp / a((f) s (x,t))ds dt, (x,£)eSM 



JO '-jo J 

is smooth. 

3. Injectivity of the X-ray transforms 

3.1. Proof of Theorem A. Now we proceed with the proof of Theorem A. Let first 
/ £ C£°(M mt ), and assume that I a f = 0. Then by Lemma [2.31 the function 

fr+(x,£) „ ft - 

u(x,£)= / /(7(t))exp / a(7^(s),7^(s))ds dt, 

JO l J0 J 

is smooth in SM and satisfies 

Gu + au = —f in SM, u\ssm = 0. 

By Theorem 12.21 there is a holomorphic w £ C°°{SM) and antiholomorphic w £ 
C°°(SM) such that 

G(e~ w u) = -e~ w f, G(e-™u) = -e~*f. 

Now Lemma I2T31 shows that e~ w u is holomorphic and e~ w u is antiholomorphic. Multi- 
plying by e w and e w , it follows that u itself is both holomorphic and antiholomorphic. 
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This is only possible when u does not depend on £, that is, u = uq. Now the transport 
equation reads 

du (£) + au = -f in SM, u \asM = 0. 
Evaluating this at ±£ and substracting the resulting expressions gives that duo = 
— u a, hence u = by the boundary condition. Consequently / = 0. 

It remains to prove the result when / G C°°(M) may have support extending up to 
the boundary. In fact, this case can be reduced to the result for compactly supported 
functions by using the general principle that (I a )*I a is an elliptic pseudodifferential 
operator. 

Suppose / G C°°(M) and I a f = 0. We consider more generally the weighted ray 
transform with weight p G C°°(SM), 

Vfo = p(Mx, ° <kix, 0)dt, (x, e <9 + sm. 

Jo 

With the choice p = e u ° we obtain 7 p / = 0. Let (M,g) CC (M,g) be another simple 
manifold which is so small that any M-geodesic starting at a point of d^SM never 
enters M again. We extend a to SM as a smooth function and / by zero to M, and 
denote by I p the corresponding weighted ray transform in M. 

It is easy to see that I p (x, £) for all (x, ^) G d + SM, since either the geodesic starting 
from (£, £) never touches M or else I p (a;, y) = I p f(x,£) for some (x,£) G d + SM. By [HI 
Proposition 2], since p is nonvanishing, I*I p is an elliptic pseudodifferential operator 
of order —1 in M int . Now I*I p f = 0, and elliptic regularity shows that / is smooth. 

Thus / G C c °°(Af int ) and I p = 0, showing that I a f = 0. The result above implies that 
/ = as required. 

3.2. Proof of Theorem B. The proof of Theorem B reduces to proving the next 
result. We say that / G C°°(SM) has degree m if /&. = for |A;| > m + 1. 

Now, word for word as in [20], we explain the identification between real- valued 
symmetric m-tensor fields and certain smooth functions on SM with degree m. Given 
such a tensor / = fi 1 -i m dx n <S> ■ ■ ■ <S> dx lm we consider the corresponding function on 
SM (henceforth referred to as the restriction) defined by 

f(x,o = fn-^e i 

Then clearly 

m 
k=—m 

where fk = f-k- Moreover if m is even (resp. odd) all the odd (resp. even) Fourier 
coefficients vanish. 

Conversely suppose that we are given a smooth real- valued function / G C°°(SM) 
with degree m. Suppose in addition that if m is even (resp. odd) then f = f + (resp. 
/ = /_), i.e. all fk and parts with k < m odd (resp. even) will vanish. Since / is 
real- valued fk = /_&. For each k > 1, the function f_ k + f k gives rise to a unique real- 
valued symmetric fc-tensor F k whose restriction to SM is precisely + f k . This can 
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be seen as follows. Recall that a smooth element fk can be identified with a section 
of T*M® k hence, its real part defines a symmetric fc-tensor. (For k — 0, fo — fo is 
obviously a real- valued 0-tensor.) 

By tensoring with the metric tensor g and symmetrizing it is possible to raise the 
degree of a symmetric tensor by two. Hence XFk : = c(Fk <8> g) will be a symmetric 
tensor of degree k + 2 whose restriction to SM is again fk + since g restricts as 
the constant function 1 to SM. Now consider the symmetric m-tensor 



8=0 

It is easy to check that the restriction of F to SM is precisely /. 
The proof of Theorem B reduces for the proof of the following 

Theorem 3.1. Let (M,g) be a simple surface, and assume thatu G C°°(SM) satisfies 
Gu + au = —ip in SM with u\d(SM) = 0. If ip G C°°(SM) has degree m > 1, then u 
has degree m — 1 . 

Proof of Theorem B. Let / and be a symmetric m-tensor and m — 1-tensor field on 
M and suppose that I a (f + 0) — 0. We write 



Then u\q(sm) — 0, and also u G C°°(SM) by Lemma 12.41 Now / + (3 has degree 
m, and u satisfies Gu + au = —(/ + /?) in SM with u\q(sm) — 0. Theorem 13.11 
implies that u has degree m — 1. We let p :— —u. As we explained above, p gives 
rise to a symmetric (m — l)-tensor still denoted by p. Since G(p) + ap = f + (3, this 
implies that d s p + ap and / + agree when restricted to SM and since they are both 
symmetric tensors of the same degree it follows that d s p + ap = f + (3. This proves 
the theorem. □ 

Theorem 13.11 is in turn an immediate consequence of the next two results. 

Lemma 3.2. Let (M,g) be a simple surface, and assume thatu G C°°(SM) satisfies 
Gu + au = —ip in SM with u\q(sm) = 0. If m > and if ip G C°°(SM) is such that 
ipk — for k < —m — 1, then Uk = for k < —m. 

Lemma 3.3. Let (M,g) be a simple surface, and assume that u G C°°(SM) satisfies 
Gu + au = —ip in SM with u\q^sm) = 0. If m > and if ip G C°°(SM) is such that 
ipk = for k > m + 1, then Uk = for k > m. 

We will only prove Lemma I3.2[ the proof of the other result being completely 
analogous. 

Proof of Lemma WM Suppose that u is a smooth solution of Gu + au = —ip in SM 
where ipk = for k < — m — 1 and u\q(sm) — 0. We choose a nonvanishing function 
v G Q m and define the 1-form 



[m/2] 

F := V A*F, 



m— 2i • 




A := -u^Gu. 
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Then vu solves the problem 

(G + a + A)(vu) = -wj) in SM, vu\ d{SM ) = 0. 

Note that vip is a holomorphic function. Next we employ a holomorphic integrating 
factor: by Theorem 12.21 there exists a holomorphic w G C°°(SM) with Gw = a + A. 
The function e w v u then satisfies 

G(e w vu) = -e w vifj in SM, e w vu\ 9{SM ) = 0. 

The right hand side e w vip is holomorphic. Now Lemma [2.31 implies that the solution 
e w vu is also holomorphic and (e w vu)o = 0. Looking at Fourier coefficients shows that 
(vu)k = for k < 0, and therefore Uk = for k < — m as required. □ 
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